Landau Zener transitions in a dissipative environment: Numerically exact results 
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We study Landau-Zener transitions in a dissipative environment by means of the numerically exact 
quasiadiabatic propagator path-integral. It allows to cover the full range of the involved parameters. 
We discover a nonmonotonic dependence of the transition probability on the sweep velocity which is 
explained in terms of a simple phenomenological model. This feature, not captured by perturbative 
approaches, results from a nontrivial competition between relaxation and the external sweep. 
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Nonadiabatic transitions at avoided level crossings play 
an essential role in many dynamical processes throughout 
physics and chemistry. They have been extensively stud- 
ied both theoretically and experimentally in, e.g., solid 
state artificial atoms 1 , 0, Q , spin flips in nanomagnets 
jij], nanocircuit QED 5,Q, the dynamics of chemical re- 
actions 0] , adiabatic quantum computation H , and in 
Bose-Einstein condensates in optical lattices [9(. 

In the pure Landau Zener (LZ) problem, the dynamics 
is restricted to two quantum states coupled by a con- 
stant tunneling matrix element Ao . A control parameter 
is swept through the avoided level crossing at a constant 
velocity v. The focus is on the final occupation probabil- 
ity of the two states. This problem was solved by various 

authors independently Gum eiei- However, in any 

physical realization, a quantum system is influenced by 
its environment leading to relaxation and phase deco- 



herence during time evolution 1J]. Any manipulation 
scheme using LZ transitions should thus include such ef- 
fects. 

Many limiting cases for the dissi pat ive LZ problem are 
2, EE E3, El, El- It has been, 



analytically tractable 
e.g., shown that for fast sweeps, at low temperatures 
or for strong system-bath couplings, the bath does not 
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i7, m 



influence the LZ probability 

temperatures but slow sweeps, the two-level system will 
be driven towards equal occupation 15, EH- The limit 
of fast quantum noise has been considered in Ref. 
However, despite these important efforts, no exact solu- 
tion for the full parameter space is available up to date. 

In this Letter, we investigate the full parameter range 
of sweep velocities, temperatures, damping strengths and 
bath cut-off frequencies by means of the quasiadiabatic 
propagator path- integral (QUAPI) (2^, 21 1. It allows to 
include nonadiabatic as well as non-Markovian effects 
yielding numerically exact results. Most importantly, 
we identify a nonmonotonic dependence of the transition 
probability on the sweep velocity which can be explained 
in terms of a simple physical picture involving a compe- 
tition between relaxation and LZ driving. Our numerical 
results also allow to address the validity of existing per- 
turbative approaches. 



Landau Zener transitions around an avoided level 
crossing of two quantum states are described by the 
Hamiltonian (H = ks = 1) 



H LZ (t) = —cr x + —a z , 



(1) 



with the tunneling matrix element Ao and the energy 
gap between the diabatic states vt, changing linearly in 
time with sweep velocity v. Here, a x ,z are Pauli matrices 
and the diabatic states are the eigenstates (| J.) and | |)) 
of o~ z . Asymptotically at times \t\ ^> Aq/v, the diabatic 
states coincide with the momentary eigenstates of Hr,z- 
The LZ problem asks for the probability Pq of the system 
to end up in the ground state at t — +oo, having started 
in the ground state at t = — oo (the corresponding one 
for the excited state is given as 1 — Pq). Its exact solution 
dates back to the year 1932 0, EH, El, El and is given 
by 



P (v, A ) = |(t (oo)| | (-co))| 2 = 1-exp - 



ttA] 
2v 



(2) 



To include environmental fluctuations on LZ transi- 
tions, we couple Ht,z to a harmonic bath (3, yielding 



At high m=HLz{t) _^Y,Mbk + bl) + ^^lb k (3) 



with the bosonic annihilation/creation operators bk/b\. 
The bath influence is captured by the spectral func- 
tion J(u>) = 2ato exp(— u/u> c ), for which we choose 
here for definitcness an Ohmic form with the cut-off 
frequency u c and the coupling strength a [141 ]. The 



Landau-Zener probability for the dissipative problem 
P = Tr { U^ 1 1 T)(T l^ool 1X1} with the time evolution 
operator Uoo = Tcxp[— i dtH{t)\ is now a function 
not only of Ao and v, but also of a,cu c and the tem- 
perature T. In the following, we use u c = 10Ao unless 
specified otherwise. We determine P by applying the nu 



21] 



merically exact and reliable QUAPI technique [2 
which is well suited for driven problems (2ll | . 

Figure [T] shows the LZ probability P versus sweep ve- 
locity v for different temperatures and for weak coupling 
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FIG. 1: The LZ probability P for various temperatures T 
is shown for a weak system-bath coupling a = 0.0016. The 
dotted lines are guides to the eye. The solid line marks the 
coherent LZ probability Pq, while the dashed line indicates 
the high-temperature limit Psd, see text. Inset: Comparison 
with perturbative approaches for T = 4Ao, see text. 



a = 0.0016. One can distinguish a regime with large ve- 
locities, v ^> Aq, from a regime with small velocities 
(adiabatic regime) v < Aq. Little influence of the bath 
is expected at low temperatures [Til, [lj| and it vanishes 
exactly at T = (T3]. This is confirmed by our numeri- 
cal results. For small v and low temperatures, T < A , 
we find P ~ 1, and thus unmodified compared to the 
pure quantum mechanical LZ result P (solid line). For 
increasing velocity, the LZ probability decreases rapidly 
and there is hardly any temperature effect in the con- 
sidered temperature range, see Fig. [T] This observation 
agrees with results by Kayanuma and Nakayama who de- 
termined the LZ probability in the limit of high temper- 
atures, Psd = 5(1 — exp(— 7rA 2 /w)) [l5| (dashed line), 
assuming dominance of phase decoherence over dissipa- 
tion. For large velocities, Psd decreases as the pure LZ 
probability, Psd ~ Plz ~ ;>7rA 2 /u, and accordingly no 
sizable temperature effect is expected. 

Beyond these limits, in the regime of intermediate to 
high temperatures, T > A and small sweep velocities, 
v < Aq, we find a nontrivial and unexpected behavior 
of P. Besides an overall decrease of P with increasing 
temperature, we find (at fixed T) for decreasing velocity 
first a maximum of P at v max < A 2 
Wmin and finally again an increase (e.g., for T 



max ,0 

laq, then a minium at 
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Vmax ~ 0.5Ao,i> m m ~ 0.02Ag). For decreasing tempera- 
tures, v m - m decreases, and P(v min ) increases. 

This nonmonotonic feature cannot be described in 
terms of perturbative approaches. Ao and Rammer de- 
rived temperature-dependent corrections to the LZ prob- 
ability for low temperatures [ljj ]. Their result is shown 
in Fig. [1] (inset, solid blue line). They also found an on- 
set temperature T D on 1/v, above which the temperature 
affects P. Thus, at larger velocities, the decrease of P 
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2: Same as Fig. [1] for larger system-bath coupling a 
(a) and a = 0.2 (b). 



due to increasing temperature starts at higher tempera- 
tures. This is in line with our findings of the maximum 
in P(v), but it does not account for the minimum and 
the subsequent increase of P for smaller v. In the limit of 
high temperatures, Psd — ► 1/2 for v < Aq. Thus, Psd 
captures the decrease of P with increasing temperature, 
but does not account for the nonmonotonic behavior for 
decreasing v. Similarly, we were not able to match nei- 
ther Eq. (54) of Ref. |l5| nor Eq. (40) of Ref. with 
our exact data, see Fig. [TJ inset. 

In the following, we provide a simple physical picture 
of how the bath influences the LZ probability. We assume 
that the bath mainly induces relaxation. Since initially 
the system is in the ground-state, only absorption can 
occur, if an excitation with energy At = yAo + (vt) 2 
exists in the bath spectrum and is thermally populated. 
Since At is (slowly) changing with time, relaxation can 
only occur during a time window \t\ < ht r (with the 
resonance time t r = 2^/ A 2 . — Aq/w), in which the en ergy 
splitting fulfills the condition A t < A c = min{T, uj c } [221 ] . 
In order for relaxation processes to contribute, the (so far 
unknown) relaxation time r r must be shorter than t r . 

For large sweep velocities, t r <r n relaxation is negli- 
gible and no influence of the bath is found as expected. In 
the opposite limit, t r 3> r r , relaxation will dominate and 
the two levels will at any time adjust their occupation 
to the momentary A t and T. Once A t > A c , relaxation 
stops since no spectral weight of the bath modes is avail- 
able and the LZ probability is P c = \[l + tanh(A c /2T)]. 
For small but finite v, equilibration is retarded, i.e., equi- 
librium is reached for an energy splitting in the past 
Af < At with t' <t and accordingly, P increases with de- 
creasing v and P(v — > 0) < P c , as observed in Fig. Q] [23| . 

Relaxation will maximally suppress the LZ transition 
when both time scales coincide, leading to a minimum of 
P at w min given by 



t r (Vmin) = T r (T,a,U c ). 



(4) 
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FIG. 3: LZ probability P for different a for T = 4A . Inset: 
Linear dependence of v m m on a. 

Within resonance, \t\ < it r , only a single phonon absop- 
tion is likely. We, thus, assume equillibration towards a 

time-averaged energy splitting A r = (2/t r ) Jq dtA t ~ 
max{A c /2, A } andP(w min ) = |[l+tanh(A r /2T)]. Sub- 
sequently, there is a maximum for a sweep velocity be- 
tween v min < w max < Aq. Fig. [5t shows P(v min ) vs. T 
which is in rough agreement with the expected behavior 
for T < Ao and T > uj Ci but the temperature indepen- 
dent range for 2A < T < uj c is not observed. 

For a fixed time and for weak coupling, we can esti- 
mate the decay rate out of the ground state using Golden 
Rule, T -1 (i) = 7ra(Ag/A t ) exp(-A t /u c )n(A t ) with the 
Bose factor n(A t ) = [exp(A t /T) - l]" 1 . For the time- 
dependent LZ problem at slow sweep velocities, we may 
assume that the bath sees a time-averaged two-level sys- 
tem and thus estimate the relaxation rate r" 1 by using 
the time-averaged energy splitting A r , i.e., 

A 2 — — 

r" 1 ~ 7ra=^ exp(- A r /w c )n(A r ) . (5) 

For increasing temperature, relaxation becomes faster, 
and, accordingly, the condition for u m ; n , Eq. ((!]), is ful- 
filled for larger velocities. This describes qualitatively the 
temperature dependence of the minimum velocity v m ; n 
in Fig. [TJ Quantitatively, this simple argument predicts 
an exponential increase at low temperatures T < 2Ao, 
a linearly decreasing v m [ n in the intermediate temper- 
ature range 2Ao ^ T < uj c , and a linear increase at 
high temperatures. Fig. \5jp shows u m i n versus temper- 
ature which nicely fits to the linear behavior u m in = 
0.019Ao+0.004A o r and thus does not show the expected 
behavior in the intermediate temperature range. How- 
ever, for T = 2Ao, we estimate v m i n = 0.028 A 2 , and for 
the linear slope in the high temperature range 0.0024Ao, 
consistent with the data. 

Clearly, this simple picture breaks down for large sweep 
velocities when the system changes its energy splitting on 
a much faster time scale than relaxation time. Thus, a 
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FIG. 4: LZ probability P for various cut-off frequencies uj c for 
T = 4Ao and a = 0.00985. Inset: P(u m in) (bottom) and » m m 
(top) vs u> c \ solid lines: predictions from our simple model. 

momentary rate and even an averaged rate is not an ap- 
propriate description anymore as it requires that at high 
temperatures and strong coupling, relaxation should be 
fast enough to partly equilibrate the system also for fast 
sweep velocities resulting in P — > i. Moreover, we notice 
that Kayanuma and Nakayama [la ] assume dominance 
of phase decoherence over relaxation in the high tem- 
perature limit. Our findings show that, whereas at very 
high temperatures or large sweep velocities, phase deco- 
herence is dominant, relaxation dominates the dynamics 
at intermediate temperatures for small sweep velocities. 

Increasing a enhances relaxation and r r decreases. 
Thus, the minimum sweep velocity is larger for a larger 
a for a fixed temperature. This is confirmed by Fig. [5] 
where P is shown for the same temperatures as in Fig. 
[TJ but for a = 0.02 (a) and a = 0.2 (b). The mini- 
mum is still observable for a = 0.02 for temperatures 
^5 T < 4Ao- At higher temperatures, only a shoulder 
remains. For a = 0.2, the local extrema disappear, but 
a monotonic growth of P with decreasing sweep veloc- 
ity is still in line with our simple picture. For coupling 
strengths a > l/v2 no bath influence is expected any- 
more [3], consistent with our data. Therefore, we focus 
on a < 0.2. We only expect a minimum (see Eq. (j4|)) for 
temperatures below the crossover temperature T c given 
by v m i n (T c ) = Aq. T c is well above the highest inves- 
tigated temperature for a = 0.0016 and well below Ao 
for a = 0.2. For the intermediate coupling a = 0.02, 
4A < T c < 12.25A and for a = 0.006 T c ~ 25A (not 
shown). T c increases roughly linearly with growing a as 
predicted by our model, but the absolute values are off 
by a factor of three. 

Surprisingly, our simple picture still holds qualitatively 
for stronger damping, when the Golden Rule is not ex- 
pected to hold. Fig. [3] shows the LZ probability for dif- 
ferent q at a fixed temperature T = 4Ao. For increas- 
ing coupling, Vmin shifts to larger velocities. In fact, 
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FIG. 5: (a) Onset temperature T and (b) maximal slope vs 
sweep velocity for a = 0.02, 0.0016. (c) LZ probability at the 
minimum sweep velocity (left axis, black circles) and mini- 
mum sweep velocity (right axis, red squares) vs temperature, 
inferred from Fig. [TJ 



«min depends linearly on a, see inset of Fig. [3] The 
linear dependence is also predicted by our model, i.e., 
Vmin = 15.8o;Aq, in very good agreement with the fit 
fmin = 17.54oAq. The decreasing maximum P(v max ) re- 
sults from the shifting minimum. Another remarkable 
fact is that the LZ probability P(v m i n ) at the minimum 
velocity is independent of a, as expected. We estimate 
the averaged splitting A r = T/2 = 2Ao, in fair agree- 
ment with A r = 2.476Ao, obtained with P(w m i n ) = 
+ tanh(A r /2T)] from the data in Fig.[3] This agree- 
ment strongly supports our conclusion that relaxation 
dominates the LZ probability in the intermediate tem- 
perature range for small sweep velocities. 

The relaxation rate (O also depends on the cut-off fre- 
quency ojc of the bath spectrum and relaxation is strongly 
suppressed when At > u> c - Fig. |4] shows P for different 
lo c , ranging down to oj c — 0.5A, a situation occurring for 
biomolecular exciton dynamics in a protein-solvent en- 
vironment [24|]. With decreasing cut-off frequencies, the 
minimum in P(v) shifts to smaller v as qualitatively ex- 
pected from Eq. ([5]) (note that this is surprising since 
a small ui c also induces strong non-Markovian effects). 
At the same time the LZ probability -P(u m m) decreases. 
With decreasing co c , the resonance time t r and the aver- 
aged energy splitting A r decrease. For cut-off frequen- 
cies T < ui c , we expect v m in = 0.19exp(— 2Aq/uj c ), in fair 
agreement with the fit i> m in = 0.23 exp(— 2Aq/lu c ). This 
is shown in the upper inset in Fig. 3J The lower inset 
shows P(v m i n ) versus w c . The solid lines are predictions 
from our model and are in fair agreement with data. 

Finally, we analyze the dependence of the onset tem- 
perature T Q on the sweep velocity v. Fig. Q] suggests that 
with increasing temperature, P decreases and the onset 
temperature increases with increasing v. We have ex- 
tracted T Q from our data by taking the intersection point 



of a linear fit of P{T) at the point of maximal slope with 
P(T = 0) = 1. Figurc[5](a) shows the onset temperature 
T and Fig. [5] (b) the corresponding slope dP/dT versus 
v. The slope dP/dT clearly displays a peak whose height 
is only weakly dependent on a and whose width strongly 
decreases with a. The onset temperature increases lin- 
early with v with a constant off-set. This finding is in 
contrast to the prediction T Q ~ 1/v of Ref. [l6| ■ 

We have investigated the dissipative LZ problem by 



means of the numerically exact QUAPI [20|, |2l| approach 
for an Ohmic bath. In the limits of large and small sweep 
velocities and low temperatures, our results coincide with 
analytical predictions [H, [l|| [l?], EH ■ For small sweep 
velocities and medium to high temperatures, however, 
we have discovered non-monotonic dependencies on the 
sweep velocity, temperature, coupling strength and cut- 
off frequency. We have shown that this parameter range 
is not accessible by perturbative means. This behavior 
can be understood in simple physical terms as a non- 
trivial competition between relaxation and LZ driving. 
This novel feature is accessible by available experimental 
techniques. 
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